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Abstract 

Directed spiral percolation (DSP), percolation under both directional and rotational constraints, 
is studied on the triangular lattice in two dimensions {2D). The results are compared with that 
of the 2D square lattice. Clusters generated in this model are generally rarefied and have chiral 
dangling ends on both the square and triangular lattices. It is found that the clusters are more 
compact and less anisotropic on the triangular lattice than on the square lattice. The elongation of 
the clusters is in a different direction than the imposed directional constraint on both the lattices. 
The values of some of the critical exponents and fractal dimension are found considerably different 
on the two lattices. The DSP model then exhibits a breakdown of universality in 2D between 
the square and triangular lattices. The values of the critical exponents obtained for the triangular 
lattice are not only different from that of the square lattice but also different form other percolation 
models. 



I. INTRODUCTION 



A new sitepercolation model, directed spiral percolation (DSP), is recently introduced 
by SantraP, |J. The DSP model is constructed imposing both directional and rotational 
constraints on the ordinary percolation (OP) model [3j. The directional constraint is in a 
fixed direction in space and the empty sites in that direction are accessible to occupation. 
Due to the rotational constraint the sites in the forward direction or in a rotational direc- 
tion, say clockwise, are accessible to occupation. The direction of the rotational constraint 
is not fixed in space and it depends on the direction from which the present site is occu- 
pied. Percolation under only directional or only rotational constraints have been studied 

n 

independently and the corresponding models are known as directed percolation (DP)[4i] and 
spiral percolation (SP)P| respectively. It is already known that both DP and SP models 
belong to different universality classes other than that of OP. The DSP model is essentially 
a combination of DP and SP models and it is constructed by imposing both the constraints 
simultaneously in the same model. Recently, the DSP model has been studied on the square 
lattice in 2 dimensions (2Z})Q,|^. It has been found that a new type of percolation cluster 
is generated in this model. They are highly rarefied, anisotropic and chiral in nature. The 
elongation of the clusters is in a different direction from the imposed directional constraint. 
The values of the critical exponents obtained are different from that of the OP, DP and SP. 
Consequently, the DSP model belongs to a new universality class. 

In this paper, the DSP model is studied on 2D triangular lattice and the results are 
compared with that of the square lattice data. The clusters are found more compact and 
less anisotropic on the triangular lattice than on the square lattice. Most interestingly, it is 
found that the values of the fractal dimension and some of the critical exponents of the DSP 
model are considerably different on the square and triangular lattices. Thus, there exists a 
breakdown of universality in the DSP model when the results of the square and triangular 
lattices are compared. In the following, the DSP model will be described briefly. The results 
of the triangular lattice will be then presented and compared with the square lattice data. 



II. THE MODEL 

Detailed description of the model on the square lattice is given in Ref . Q] . A brief demon- 
stration will be given here on the triangular lattice. A left to right directional constraint and 
a clockwise rotational constraint are imposed on the system defined on a triangular lattice 
of size L X L. Due to the directional constraint an empty site on the right of an occupied 
site and due to the rotational constraint the empty sites in the forward direction or in the 
clockwise direction can be occupied. To generate clusters under these two constraints a 
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single cluster growth algorithm is developed in Ref. jj| following the original algorithm of 
Leathjfi]. In this algorithm, the central site of the lattice is occupied with unit probability. 
All six nearest neighbours of the central site on the triangular lattice can be occupied with 
equal probability p in the first time step. As soon as a site is occupied, the direction from 
which it is occupied is assigned to it. Selection of empty nearest neighbour in the next MC 
time steps is illustrated in FigH Two long arrows from left to right in Fig^ represent the 
directional constraint. The presence of the rotational constraint is shown by the encircled 
dots. The black circles represent the occupied sites and the open circles represent the empty 
sites. The direction from which the central site is occupied is represented by a short thick 
arrow. Now the nearest neighbours of the central occupied site eligible for occupation will 
be identified. The dotted arrow indicates the eligible empty site for occupation due to the 
directional constraint and the thin arrows indicate the eligible empty sites for occupation due 
to the rotational constraint. Since the directional constraint is to the right, site 4 is always 
eligible for occupation. The rotational constraint acts in the forward or in the clockwise 
direction with respect to the direction of approach to the present occupied site. Since the 
central site is approached from 2, there are three sites, site 5 in the forward direction and 
sites 6 and 1 in the rotational direction, eligible for occupation due to rotational constraint. 
Note that, on the square lattice only two empty sites due to the rotational constraint were 
eligible for occupation at any MC step. It is also important to note that, in this model an 
occupied site can be .eoccupied fro. : differ, d.e.io,. due .o .ota.io„a, co,.Ua.tQ. 
A site is forbidden for occupation from the same direction. On the triangular lattice, a site 
then could be occupied at most 6 times from 6 different directions. Due to the reoccupa- 
tion of occupied sites, cluster generation is time consuming in the rotationally constrained 
mo dels 



After selecting the eligible sites for occupation, they are occupied with probability p. The 
coordinate of an occupied site in a cluster is denoted by {x,y). Periodic boundary condi- 
tions are applied in both directions and the coordinates of the occupied sites are adjusted 
accordingly whenever the boundary is crossed. At each time step the span of the cluster 
in the x and y directions = Xmax — and Ly = ymax — Vmin are determined. If 
or Ly > L, the system size, then the cluster is considered to be a spanning cluster. The 
critical percolation probability pc is defined as below which there is no spanning cluster and 
at p = pc a. spanning cluster appears for the first time in the system. 



III. RESULTS AND DISCUSSIONS 



Simulations are performed on the triangular lattice of several different lattice sizes from 
L = 128 to L = 2048. The cluster size distribution Ps{p) is defined as Ps = Ng/Ntot 
where is the number of s-sited finite clusters in a total of Nt^t clusters generated. The 
percolation threshold Pc at which a spanning cluster appears for the first time in the system 
is determined by generating Ntot = 5 x 10^ clusters at different site occupation probability p. 
The probability to have a spanning cluster at a given site occupation probability p is given 
by Pgp = risp/Ntot = 1 — Z^s Ps{p), where Ugp is the number of spanning clusters out of Ntot 
clusters. The percolation threshold pc is determined from the maximum slope {dPsp/dp)max 
of the curve Pgp versus p. In Fig|21 Pgp and dPsp/dp are plotted against p ioi L = 2048. 
The percolation threshold Pc could be identified as Pc = 0.5700 ± 0.0005 corresponding to 
the maximum slope. The derivative is calculated using the central difference method for the 
data points collected in an interval of 0.0005. 

An infinite cluster generated on the triangular lattice of size L = 256 at pc = 0.5700 
is shown in FiglHl The black dots are the occupied sites and the solid black circle on the 
upper left corner is the origin of the cluster. It could be seen that the elongation of the 
spanning cluster is almost along the left upper to the right lower diagonal of the lattice 
as it was seen on the square lattice. In the case of charged particles, this is due to the 
development of Hall voltage across the sample and perpendicular to the applied in-plane 
electric field. As a result, an effective directional field is developed along the left upper to 
the right lower diagonal of the lattice. It has already been observed on the square lattice 
that the DSP clusters are not merely the DP clusters in the presence of the effective field. 



The DSP clusters contain features other than the DP clusters. The clusters here contain 
holes of almost all possible sizes and it has clockwisely rotated (chiral) dangling ends. It 
could also be noticed that the infinite cluster on the triangular lattice is more compact and 
less anisotropic in comparison to the infinite cluster on the square lattice. 

The fractal dimension df of the infinite clusters at Pc = 0.5700 on the triangular lattice 
of size L = 2048 is determined by the box counting method. The number of boxes Nb{€) 
is expected to grow with the box size e as Nb{^) ~ e'^^ where df is the fractal dimension. 
In Fig. m Nsie) is plotted against the box size e. The data are averaged over 5 x 10^ 
spanning clusters. A reasonably good straight line is obtained in the log — log scale. The 
fractal dimension is found df = 1.775 ± 0.004. The error is due to the least square fitting of 
the data points taking into account the statistical error of each point. In order to check the 
convergence of the value of the fractal dimension, df is plotted against l/N^p in the inset of 
Fig. H It could be seen that the value oi df remains unchanged over 10^ to 5 x 10^ spanning 
clusters. The value of df has also been estimated from finite size (FS) scaling So^ ~ L'^^ , 
where 5*00 is the size of the largest cluster aXp = pc- The lattice size L changes from 2^ to 2^^. 
It is found that df{FS) = 1.80 ± 0.03, which is within the error bar of the other estimate. 
The fractal dimension df obtained here is higher than that of df ^ 1.733[lJ {df{FS) ^1.72 
[2 1) on the square lattice. Also notice that the value of df obtained here is smaller than the 
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fractal dimensions obtained in OP (91/48 8]) and SP (1.969 and it is slightly higher 
than DP (1.765 [7J). Vacancies are generated into the cluster as it grows. At the same time, 
due to the higher number of branching on the triangular lattice the cluster penetrate into 
itself more and more than on the square lattice. As a result, the infinite clusters are less 
rarefied on the triangular lattice than on the square lattice. 

Since the fractal dimension df is different from that of the square lattice value, it is then 
expected that the values of the other critical exponents will also be different form that of 
the square lattice in order to satisfy the scaling relations among the critical exponents. The 
critical exponents related to the different moments of the cluster size distribution Ps{p) are 
now estimated. The scaling function form of the cluster size distribution Ps{p) for single 
cluster growth technique, in which the central site is occupied with unit probability, is 
assumed to be 

P,(p) = ,-+if[s-(p-p,)] (1) 
where r and a are two exponents. The assumed form of the scaling function Ps{p) is the same 



as that of the square lattice. The first moment x = Yl!s^Ps{p) corresponds to the average 
cluster size. Next two higher moments are defined as xi = YHs s^Psip) and X2 = J2's s^Ps{p)- 
The primed sum represents the sum over all the finite clusters. As p ^ pc, the moments x, 
Xi, and X2 of Ps{p) become singular with their respective critical exponent 7, 6, and rj. The 
critical exponents are defined as 

x^b-Pcl"^, xi^\p-Pc\~\ & X2 ~ b-j^d"''- (2) 

To estimate the values of 7, 5 and r] on the triangular lattice, the average cluster size x and 
two other higher moments Xi and X2 are measured generating 5 x 10^ finite clusters below 
Pc for different p values on several lattice sizes. In FiglHl X; Xi and X2 are plotted against 
\p — Pc\ for the system size L = 2048. The circles represent X; fhe squares represent xi 
and the triangles represent X2- The values of the exponents obtained are 7 = 1.98 ± 0.01, 
S = 4.30 ± 0.02 and 77 = 6.66 ± 0.04 for L = 2048. The errors quoted here are the standard 
least square fit error taking into account the statistical error of each single data point. 
Because of the error bar Ap^ = 0.0005 in the threshold, all the exponents have also been 
estimated for two other critical probabilities pc^Apc- The values of the exponents obtained 
for p = 0.5695 are 7 ~ 1.96, 6 ~ 4.27, and rj ^ 6.60 whereas for p = 0.5705 they are 
7 ^ 2.00, S ~ 4.33 and rj ^ 6.71. The values of the critical exponents are then taken as: 
7 = 1.98 ± 0.02, 6 = 4.30 ± 0.04 and r] = 6.66 ± 0.08. The values of the critical exponents 
at the optional thresholds Pc ± Apc are now within error bars. A comparison of the values 
of the critical exponents obtained on the triangular and square lattices is made in FigEl 
for several lattice sizes. In FiglHl the values of 7, 6 and 1] are plotted against the inverse 
system size 1/L. The squares represent the square lattice data and the triangles represent 
the triangular lattice data. The data for the square lattice is taken from Ref.P] except 
for L = 2048. In Ref.[l|, data were reported upto the maximum lattice size L = 1024. 
For the sake of comparison with L = 2048 triangular lattice data, new estimates of the 
critical exponents have also been made on L = 2048 square lattice. The results of L = 2048 
square lattice are in good agreement with that of the smaller system sizes. On the square 
lattice of size L = 2048, the values of the critical exponent obtained are: 7 = 1.85 ± 0.01, 
S = 4.01 ± 0.04 and 77 = 6.21 ± 0.08. The results on both the lattices are then extrapolated 
upto L 00, the infinite system size. The extrapolated values of the exponents are marked 
by crosses. The values of the critical exponents are found very different (beyond the error 



bars) on the square and triangular lattices. The triangular lattice values of the exponents 
are higher than that of the square lattice. The clusters then grow much larger in size on 
the triangular lattice than on the square lattice for a given p. This might be due to higher 
number of branching possibilities on the triangular lattice. Also notice that the value of 
2(5 — 7 = 6.62 is very close to the value of the exponent rj = 6.66. The values of the 
exponents then satisfy the scaling relation rj = 26 — 'j [1] within error bars. The exponents 
are not only different from the square lattice values but also different from that of other 
percolation models, OP Q , DP Q and SP Q • 

There are two connectivity lengths, and for the anisotropic clusters. Here, 
is along the elongation of the cluster and is along the perpendicular direction to the 
elongation. The connectivity lengths are defined as = 2^2'^ R'jsPs{p)/ sPs{p) and 

= 2^2'^ R^sPsi^p) / sPs{p) where and R± are the radii of gyration with respect to 
two principal axes of the cluster. They are estimated from the eigenvalues of the moment of 
inertia tensor, a 2 x 2 matrix here. C,\\ and diverge with two different critical exponents 
z/j| and z/_L as p — s> Pc- The critical exponents z/y and are defined as 

^ii-b-Pcl""" & a-b-Pcl""^. (3) 

The connectivity lengths, ^y and for the system size L = 2048 are plotted against \p~Pc\ 
in Fig|7| Data are averaged over 5 x 10^ clusters. The squares represent ^y and the circles 
represent C,±- The corresponding exponents are z/y = 1.36 ± 0.02 and z/^ = 1.23 ± 0.02. 
The errors quoted here are the least square fit errors. The values of the exponents are also 
estimated at pc ± Apc- For p = 0.5695, the values obtained are z/y ~ 1.35 and z/_l ^ 1.21 
and for p = 0.5705, the values obtained are z/y ^ 1.37 and i>± ^ 1.24. There is a little 
variation and the values of the critical exponents are approximated as: z/y = 1.36 ± 0.02 
and z/x = 1.23 ± 0.02. The error bars now include the values of z/y and z/^ at the optional 
PcS. To compare the square and triangular lattice data, simulations have been performed 
on other smaller system sizes. In Fig. |H1 the exponents z/y and z/_|_ are plotted against the 
inverse system sizes 1/L for both the square and triangular lattices. The squares represent 
the square lattice data and the triangles represent the triangular lattice data. Data of the 
square lattice is taken from Ref. Q except for L = 2048. New estimates of the exponents are 
also made on L = 2048 square lattice. The exponents are extrapolated upto L oo and 
they are marked by crosses. Notice that, the exponent z/y is almost the same as that of the 



square lattice value (~ 1.33) whereas is higher than that of the square lattice (~ 1.12). 
Both the exponents are also different form that of DP model l|, llO]. The hyperscaling 
relations 25 — 37 = (d — 1)z/_l + z/y and {d — df)u± = /5 = 5 — 27 jl|] are satisfied marginally: 
25 - 37 = 2.66 ± 0.06 whereas {d - l)z/^ + u\\ = 2.59 ± 0.04 and 5 - 27 = 0.34 ±0.06 whereas 
{d — df)i>± = 0.28 ± 0.03. It is already known that the hyperscaling is violated in directed 
percolation ^|. In the case of DSP, the hyperscaling relations were found satisfied on the 
square lattice whereas on the triangular lattice they are "marginally" satisfied. The ratio of 
the connectivity lengths goes as ~ \p — Pc\~^'^ where Au = z/y — z/_l. For the square 

lattice, Az/ is approximately 0.21 whereas for the triangular lattice, it is approximately 0.13. 
Thus, the clusters are less anisotropic on the triangular lattice. This is because of more 
flexibility given to the spiraling constraint which makes the cluster not only compact but 
also less anisotropic. 

The values of the critical exponents and fractal dimension obtained in the above study for 
the triangular lattice are summarized and compared with the square lattice data in TableHl 
The values in the parenthesis are the suggested rational fractions for the values of the critical 
exponents on the square lattice. These rational fractions satisfy the scaling relations exactly 
including the hyperscaling relations. It can be seen that some of the critical exponents and 
the fractal dimensions are considerably different on the square and triangular lattices for 
the DSP model. According to the theory of critical phenomena, the values of the critical 
exponents are independent of the underlaying lattice structure in the same spatial dimension. 
As a consequence, the systems defined on different lattices in the same space dimension then 
belong to the same universality class. Since the values of the critical exponents of the DSP 
model differ on the square and triangular lattices in 2D, the DSP model then exhibits a 
breakdown of universality. This is the first percolation model which shows breakdown of 
universality on two different lattices in the same spatial dimension. It is already seen in the 
above discussion that the flexibility in the spiraling constraint makes the clusters compact 
and less anisotropic. A possible reason for different critical behaviour on the square and 
triangular lattices may be due to different scaling behaviour of the finite clusters below 
percolation threshold on the two lattices. Below Pc, the finite clusters are called lattice 
animals jl^. Lattice animals without any loop are known as lattice trees. Though the 
spiral lattice animals have the same scaling form on the square and triangular lattices, 
it has been found that the spiral trees (lattice animals without loops) follow two different 



scaling relations on the square and triangular lattices and belong to two different universality 
classes^]. In the asymptotic n ^ oo limit, the number of spiral lattice site trees (a„) of 
n-s,tes on the triangular latt.ce obey the scahng plation g,ven by a„ « A"'.- Q whereas 
on the square lattice it is given by a„ ^ A^ra ^ where 6 and 6 are two exponents and 
A is known as the growth parameter. The origin of different scaling forms for the spiral 
lattice site trees on the square and triangular lattices is due to the fact that they can not 
have branching on the triangular lattice except at the origin whereas on the square lattice 
branching is possible at any point. The radius of gyration exponent of spiral trees was also 
found different on the square ~ 0.653) and triangular lattices (z/ ^ 0.618) Another 
lattice statistical model, the spiral self-avoiding walks (SAW) also exhibit breakdown of 
universality on the square and triangular lattices. The asymptotic (large n) behaviour of 
the number of walks Sn is given by Sn ~ An~'^ exp{Xy/n) , where A = 2' -2 X 3-5/^7r, 7 = 7/4, 
and A = 2^l/^/3 for the square lattice Q and A = 2^/^ x 3-^/^7r, 7 = 5/4, and A = n/,j2/3 
for the triangular lattice jl^. Notice that the scaling relation for the spiral lattice site 
trees is similar to that of the spiral SAWs. It should be mentioned here that the values of 
the critical exponents and the scaling behaviour of the cluster related quantities in the SP 
model (percolation in the presence of rotational constraint only) are the same on the square 
and triangular lattices and on breakdown of universality has been observed This may 
be due to the fact that the spiral lattice trees are minority in number at the percolation 
threshold and unable to change the universality class. In the DSP model, the presence of the 
directional constraint on top of the rotational constraint might increase the number of spiral 
lattice trees. Maybe, the higher number of tree like structures in the clusters generated has 
a non-trivial effect on the critical properties of the DSP model at the percolation threshold 
and leads to breakdown of universality of the critical exponents. 

Finally, the form of the scaling function Ps{p) = s~'^+^f [s°"(p — pc)] is verified. The 
exponents r and a are estimated using the scaling relations /3 = (r — 2)/cr, 7 = (3 — r)/cr, 
6 = (4 — r)/cr, and rj = (5 — r)/cr and following the same technique described in Ref. . The 
estimates of r and a are obtained as r = 2.16 ±0.02 and a = 0.427 ±0.003 respectively. The 
errors quoted here are the propagation errors. On the square lattice, the values of r and 
a were obtained as ~ 2.16 and 0.459 respectively. Note that, the value of r is the same as 
that of the square lattice whereas a differs on the two lattices. The scaling function form is 
verified through data collapse by plotting Ps{p) / Ps{pc) against the scaled variable s'^{p — p^ 



in Fig. 121 The cluster size s changes from 64 to 16384 and {p — pc) varies from 0.007 to 
—0.06. A reasonable data collapse is observed. The scaling function form is similar to that 
of the square lattice. The height of the function remains almost the same but the width 
is slightly smaller than that of the square lattice function. The lesser width of the scaling 
function Ps{p) is just a consequence of the lesser value of a on the triangular lattice. 

IV. CONCLUSION 

The directed spiral site percolation is studied on the triangular lattice and the results 
are compared with that of the square lattice. Clusters on the triangular lattice are found 
more compact and less anisotropic than the clusters on the square lattice. Interestingly, it 
is also found that the values of the fractal dimension and some of the critical exponents on 
the triangular lattice are significantly different from the square lattice values. This might be 
due to different scaling behaviour of some of the finite clusters below percolation threshold 
on the two lattices. As a consequence, the DSP model exhibits a breakdown of universality 
between the square and triangular lattices in 2D. The values of the critical exponents on 
both the lattices satisfy the scaling relations between the moment exponents (7,5,77). The 
hyperscaling relations were satisfied on the square lattice but they are "marginally" satisfied 
on the triangular lattice. The exponents are not only different on the square and triangular 
lattices but also different from other percolation models like OP, DP and SP. Directed spiral 
percolation is expected to occur in disordered systems when both rotational and directional 
force fields are present. 
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1.72 ± 0.02 (FS) 

Triangular: 1.775 ± 0.004 1.98 ± 0.02 4.30 ± 0.04 6.66 ± 0.08 1.36 ±0.02 1.23 ±0.02 
1.80 ±0.03 (FS) 

TABLE I: Comparison of the critical exponents and fractal dimension of the DSP model measured 
on the square and triangular lattices. For the square lattice, the values within parenthesis are the 
suggested rational fractions for the values of the critical exponents in Ref.Q| which satisfy the 
scaling relations exactly. Some of the critical exponents and the fractal dimension are significantly 
different on the two lattices. 




FIG. 1: Selection of eligible nearest neighbors for occupation of an already occupied site. Black 
circles are the occupied sites and open circles are the empty sites. Two thick long arrows from 
left to right represent the directional constraint. The presence of clockwise rotational constraint 
is shown by the encircled dots. The eligible nearest neighbors of the central occupied site will be 
selected here for occupation. Six nearest neighbours of the central site on the triangular lattice 
are marked as 1 to 6. The central site is occupied from the site 2, marked by a thick arrow. Due 
to directional constraint, site 4 on the right of the occupied site, is always eligible for occupation 
and it is indicated by a dotted arrow. Due to rotational constraint, sites 5, 6, and 1 are eligible for 
occupation and they are indicated by thin arrows. 
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FIG. 2: Plot of spanning probability Pgp and the slope dPsp/dp versus p. The circles represent 
Pgp and the squares represent the slope dPgp/dp. The critical probability Pc is determined from 
the maximum slope. For the triangular lattice, it is found that Pc = 0.5700 ± 0.0005 as indicated 
by an arrow. 



FIG. 3: An infinite cluster on a 256 x 256 triangular lattice at p = 0.5700 is shown. The black dots 
are the occupied sites. The solid black circle on the upper left corner is the origin of the cluster. 
The cluster has holes of almost all possible sizes. The elongation of the cluster is along the upper 
left to the lower right diagonal and not along the imposed directional constraint from left to right. 
The dangling ends are clockwisely rotated. 
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FIG. 4: Number of boxes Nsie) is plotted against the box size e. Data are averaged over 50, 000 
spanning clusters generated on a triangular lattice of size L = 2048. The fractal dimension is found 
df = 1.775 lb 0.004. In the inset, df is plotted against l/Ngp, the number of spanning clusters. It 

could be seen that the value has converged with the realizations. 
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FIG. 5: Plot of the first, second and third moments Xi Xi) ^-nd X2 of the cluster size distribution 
versus |p — Pc| for a triangular lattice of size L = 2048. Different symbols are: circles for x, squares 
for xi, and triangles for X2- The solid lines represent the best fitted straight lines through the data 
points. The corresponding critical exponents are found as 7 = 1.98 it 0.02, 5 = 4.30 it 0.04, and 
r? = 6.66 ± 0.08. 
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FIG. 6: Plot of the exponents 7, 5 and rj against the inverse system size 1/L. The system 
sizes considered are: L = 128, 256, 512, 1024 and 2048. The squares represent the square lattice 
data and the triangles represent the triangular lattice data. The exponents are extrapolated upto 
L ^ 00 and the extrapolated values are marked by crosses. It can be seen that the exponents are 
significantly different on the two lattices. 
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FIG. 7: The connectivity lengths, .^|| and are plotted against |p — ^icl for a triangular lattice of 
size L = 2048. The circles represent and the squares represent ^y. The solid lines represent the 
best fitted lines through the data points. The critical exponents are found as i/y = 1.36 it 0.02 and 
= 1.23 ± 0.02. 
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FIG. 8: The connectivity exponents and arc plotted against the inverse system size 1/L for 
the square and triangular lattices. The system size changes form L = 128 to 2048 as in Fig.6. The 
squares represent the square lattice data and the triangles represent the triangular lattice data. 
Extrapolated values to the infinite system size (1/L = 0) are marked by crosses. The value of 
seems to be different on the two lattices whereas i^ii is close to the square lattice value. 




FIG. 9: Plot of the scaled cluster size distribution Ps{p) / Ps{pc) versus the scaled variable s'^{p—pc) 
for different values of p on the triangular lattice. The value of a is taken as cr = 0.427. The cluster 
size s changes from 64 to 16384. The data plotted correspond to p — Pc = 0.007(x), 0.005(n), 
-0.035(v), -0.04(0), -0.045(<), -0.05(*), -0.055(A), -0.06(0)- A reasonable data collapse is 
observed. 



